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Abstract. This review concentrates on the results obtained, over the last ten years, on the 
astrophysics of high-energy cosmic ray electrons and positrons. The anomalies, observed in 
the data of recent experiments (possible bump in the electron spectrum and the PAMELA 
anomaly in the positron fraction) are discussed through the systematic use of simple analytical 
solutions of the transport equations for cosmic ray electrons. Three main ways of explaining the 
origin of the anomalies are considered: the conservative way supposing the positrons to be pure 
secondary particles; the nearby sources like pulsars origin; and the dark matter origin. This 
review discusses, also, the inability to select the pulsars model or the dark matter model to 
explain the electron anomalies on the basis of the electron spectra with the usual large energy 
binning (> 15%). It is argued that the signature of nearby pulsars origin of the anomalies 
against the dark matter origin could be the fine structure of the cosmic ray electron spectrum 
predicted in the Malyshev et al. paper (2009) and which was observed in the data from the 
high-resolution ATIC experiment (2009-2011). To date, the high-resolution ATIC data was the 
only experimental result of this type published in the literature. Therefore, they should be 
tested by other experiments as soon as possible. Generally, there is, also, rather controversial 
situations between the data of the majority of recent experiments and, consequently, there is a 
noted urgent need for new high-precision and high-statistical experiments. 



1. Introduction 

Over the last few years a large amount of new high-precision experimental data related to primary 
electrons and positrons in cosmic rays was obtained in the experiments of a new generation: 
PPB-BETS P, ATIC 0131, PAMELA gUSEl, Fermi/LAT [2 El [9], H.E.S.S. HQ], MAGIC p]. 
Some observed features of the data were unexpected and, now compared to the beginning of the 
2000's, the situation looks very exciting and intriguing. The difference can be seen clearly if one 
compares the content of the present review with D. Mtiller paper |12] with an almost similar 
titkQbut published in 2001. The main subject of present review is the total energy spectrum of 
cosmic ray electrons and positrons in the energy range higher than approximately 100 GeV, the 
positron fraction in the total flux, and related astrophysics. A number of interesting problems 
like charge-dependent solar modulation; anisotropy of electron and positron fluxes; and some 
others are not considered. Hereinafter, in this paper, the term "electrons" is used for the sum 
of all cosmic ray charged leptons regardless of the charge. Otherwise, "negative electrons" or 
"positrons" are used. 

In the total cosmic rays flux, the fraction of electrons is small (see Figure [1]). However, due to 
their special properties, they are very important to the astrophysics. For high-energy electrons 
{E > 10 GeV), the dominant channel of energy losses is synchrotron radiation due to interaction 



1 "Cosmic-ray electrons and positrons" . 
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Figure 1. Ratio of the electron flux to the 
flux of protons given by the ATIC data. The 
ratio was obtained with the proton spectrum 
from the paper [13] and with the ATIC 
electron spectrum shown in Figure El the 
present paper. 
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Figure 2. To the thin disk approximation - 
the geometry of the Galactic disk and diffusion 
halo. 



of electrons with interstellar magnetic fields and inverse Compton scattering. In the Thompson 
approximation, that holds for electrons very well up to energies around a few TeV, the rate of 
the energy loss is given by: 

^ = -b E* EE b(E) (1) 



where the factor 6q is 
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In ([2]) W crn b, Wdust, W s tar are the energy density of the cosmic microwave background (CMB), 
infrared, and starlight photons respectively, B is the mean-square induction of the magnetic 
field, and ax is the Thompson cross section. Due to the small mass of electrons compared to 
the proton mass, the radiative energy losses of electrons are much stronger than the radiative 
losses of protons and nuclei. The value of bo for electrons is not known exactly but it is believed 
widely that b = (1.2 -r 1.6) • 10" 16 (GeV • s) _1 . 

The solution of the equation (TJQ) with the initial condition E(0) = Eq is 

E(t) = ^r=-. (3) 

It can be seen easily from ([3]) that the electrons with the energy Eq could not be older than 

(E ) = l/(b E ) (4) 

since E{t) — > oo as t — > —1/(PqEq). For Eq = 1 TeV the formula (HJ) produces the value 
Tmax ~ 2 • 10 5 years which is much shorter than the lifetime of protons in the Galaxy 
(~ 10 7 years). In the diffusion approximation the mean distance which could be passed by 
high-energy electrons during their lifetime is determined approximately by the diffusion length 



Tmax ~ V ^DT ma x • 



(5) 



For the diffusion coefficient D, the widely adopted expression is: 

D = D ■ (E/GeY) s , (6) 

where 5 = 0.3 -r 0.6 and D = (1 -r 5) • 10 28 smV 1 . Then for b = 1.4 • 10~ 16 (GeV • s)" 1 and 
with formula © , eq. (0) may be written as 



6.4 / 
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For D = 3 • 10 28 smV 1 , one obtains from flTJ): r moa; (l TeV, 5 = 0.6) = 1.5 kpc, r max (l TeV, 5 = 
0.3) = 0.6 kpc. It can be seen that the flux of the ~TeV-energy electrons are sensitive to 
nearby sources in the range of the order of only one kiloparsec and this is very different from the 
nuclei component of cosmic rays. It was recognized clearly more than forty years ago [14j . The 
sensitivity of electron flux to a few nearby sources like pulsars or supernova remnants (SNR) 
makes it very important for astrophysics. The other reason of the importance of cosmic ray 
electrons for astrophysics is a possible connection of them with dark mater (DM) annihilation 
or decay and with other exotic sources like black hole evaporation etc. 

There are numerous more or less conservative or, on the contrary, very exotic scenarios and 
models to describe and understand the data related to the cosmic ray electrons. It is convenient 
to start with the most conservative one - the so called 'conventional model'. 



2. Diffusion approximation, thin disk approximation, and conventional model 

2.1. The diffusion equation and Green function solution 

Since S.I. Syrovatskii seminal paper [15], the diffusion approximation is considered to be a 
reasonable approximation in the transportation problem of cosmic ray electrons. In the simplest 
case of homogenious and isotropic diffusion, the diffusion equation reads: 



dp(r,E,t) d 



[b(E)p]+D(E)V 2 p + Q(r,E,t). 



(8) 



dt dE 

Here, p(r,E,t) is the density of electrons of energy E; b(E) defined by eq. ([JJ; and Q(r,E,t) is 
the source function of electrons. Equation ([S]) is linear and, therefore, for an arbitrary source 
function, the solution may be written with Green function G(r, t, E\r' , if, E'): 

p(r,t,E) = J G(r,t,E\r',t',E')Q(r',t',E')d 3 r'dt'dE', (9) 

where the Green function is [15J: 
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is the time of cooling of an electron from energy E' to E and 
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has the meaning of the mean square of diffusion length during the cooling of an electron from 
E' to E. Hereinafter in this paper we measure E in GeV and do not write explicitly (E/GeV) 
in formulae. 



2.2. Solution of diffusion equation in thin disk approximation 

In this review we use various simplified analytical solutions of the transport equation ([8]) 
systematically to understand and illustrate the physics of cosmic ray electrons. Whilst such 
models are not exact, they can be very useful and instructive. One such model is the model of 
infinitely thin and homogeneous Galactic disk together with an infinitely thick Galactic halo. 

The half-depth of the diffusion Galatic halo is ~ 4 kpc which, for 100 GeV-1 TeV electrons, 
is larger than the expected diffusion electron range r max ~ 1 kpc. At the same time, the half- 
depth of the Galaxy disk, at the position of the Sun, is only about 150 pc which is much smaller 
than r max (see Figure [5]). Therefore in first approximation, for the sources of electrons located 
within the Galatic disk, the depth of the halo may be considered as infinitely large. At the same 
time, the source of electrons located within the Galatic disk, may be considered to be infinitely 
thin relative to the value of r max . Moreover, the Sun is located very close to the Galatic plane, 
therefore, in this model, the flux of electrons, calculated exactly for z = 0, is a reasonable 
approximation in understanding the electron spectrum in the energy range 100 GeV-1 TeV. 

Now, we consider a source of electrons with spectrum Q(E), constant in the time, and 
distributed homogeneously at z = in the infinitely thin plane. Then, the source function 
reads: 

Q{r,t,E) = Q(E)5(z) (13) 
and the solution ([9]) of the transport equation ([8]) for z = reduces to: 

For power-law spectrum Q(E) oc E' 1 and with diffusion coefficient, defined by ([6]), the integral 
in eq. (fl~4"|) is calculated easily and, for the observed spectrum, one obtains: 

p{E)\ z=Q = Q,E~^- A = 5+i, (15) 

where Qo is some constant (could be calculated explicitly). Therefore, instead of the source 
spectrum with the index 7, an observer measures the electron spectrum steeper on A = S + 1/2. 
Since 0.3 < 5 < 0.6, then 0.8 < A < 1.1. In other words, Awl. It is a rather robust result. 
It does not depend on the model fine details. For example, for a model of infinite homogeneous 
source, instead of a thin plane source, one checks easily that the solution for power-law spectrum 
would be p{E) = Q E~^ +1 \ It is rather close to {15]) since A ~ 1. Please note that, throughout 
the paper, we neglect the effects of Solar modulation since these may be important only for the 
energies below ~ 10 GeV, however we are interested mainly in energies above 10-100 GeV. 

2.3. "Conventional model" in the thin disk approximation 

Cosmic ray electron spectrum, which was calculated in I.V. Moskalenko and A.W. Strong paper 
|16j , are referred commonly to as "conventional model'H. A family of similar models were studied 
also earlier by R. J. Protheroe in [TP]. The conventional model is the most conservative approach 
in understanding the cosmic negative electron and positron fluxes. In the model, it is supposed 
that there exist two main sources of negative electrons and positrons. 

Generally, the first kind of sources is the same for primary negative electrons and nuclei 
(like supernova remnants, SNR). However, sources of this kind do not accelerate positrons. 
In the conventional model, this kind of sources is considered in continuous and homogeneous 

2 Also, several variants of similar calculations with various suppositions were published by the same authors later 

HHEOU- 



(14) 



approximation. The primary negative electrons source spectral index and the intensity of the 
source were tuned in [16] to fit the experimental dated- 

The second kind of source is secondary electrons produced in inelastic scattering of cosmic 
ray protons and other nuclei in interstellar medium (ISM) gas. Moreover, it is supposed in the 
conventional model that this kind of the source is the only source for positrons. 

The spectra were calculated in [16] by numerical solution of diffusion transport equation for 
a realistic distribution of the substance in the Galaxy. This was done by using the well-known 
GALPROP system [22] developed by the authors of paper [T5]. However it is meaningful to 
study the conventional model by a more simple analytical solution of the transport equation 
such as homogeneous thin disk approximation, eqs. (|14U15|) since such models show very clearly 
the underlying physics and main features of the result. 

We suppose power-law source function with spectral index 70 for primary negative electrons 
and no positrons at all for SNR as a source. The spectral index 70 is not known exactly but it 
is commonly believed 2.1 < 70 < 2.5 and it is expected 70 to be close to the spectral index of 
proton source spectrum since the sources of the primary electrons and cosmic ray protons are 
generally the same [23\ 112] . The index 70 should be considered as a matter of fit of a model 
to the data. The main source of the secondary electrons is cosmic ray protons that produce 
them via ir^ decay after interaction with ISM nuclei. The spectrum of secondary electrons is a 
power-law spectrum with the spectral index equal to the spectral index of the observed spectrum 
of cosmic ray protons, 7 S ~ 2.7. Then eq. (|15p implies the total spectrum of negative electrons 
together with positrons to be: 

Qe- + e + (E) = Q ^- (70+A) + Q s ,e-£- (7s+A) + Qs,e + E~^ +A) , (16) 

where Qo, Q s e -, Q s ^ e + are amplitudes of the primary negative electron flux, secondary negative 
electron flux, and secondary positron flux respectively. The fit of the conventional model to the 
data in the paper [16] shows that the flux of the primary electrons is much higher than the flux 
of the secondary ones. Therefore the first term in the sum ([16]) dominates strongly and, actually, 
the shape of the spectrum Q e - +e +(E) is close to the power law with the spectral index (70 + A). 
For the fraction of positrons in the total electron flux one obtains easily: 



Qe+{E) = Qs 

Q e - +e+ (E) Q s>e - + Q Si6+ + Q £(7s-7o) • 



R e+/(e~+e+){ E ) = 7^ TFa = ~ ~ ^ W^-^V ( 17 ) 



The value (7s — 7o) is essentially positive in its nature since it is supposed that primary electrons 
and protons have similar source spectrum indexes and 7 S ~ 70 + 5. Therefore, the ratio 
Re+/(e-+e+)(E) ls a decreasing function of energy. 

Thus, an approximate power-law spectrum with the spectral index (70 + A) > 3 for the total 
electron flux, and a decreasing fraction of positrons are generic predictions of the conventional 
model. These features do not depend on details of the conventional model realization. 

The conventional model is used commonly as a reference point for discussion of electron data. 
Also, from the point of view of the conventional model being the first step of the analysis, we 
present and discuss the data from the experiments. 



3. Electron spectrum and positron fraction: the current state of the experimental 
studies 

3.1. The electron spectrum 

In the total electron spectrum, the most interesting features, measured in recent experiments, 
were observed above 100 GeV. Consequently, we restricted ourselves to reviewing only the data 



3 The most important sources of the data for paper [TS] were papers [20] and [21] . 




Figure 3. Measuremets of (e + + e~)- 
spectrum before 2000 and the conventional 
model. The experiments, presented in the 
figure are: Zatsepin, 1971 [M]; Meyer, 1971 
[27]; Silberberg, 1973 [M]; TIFR(Anand), 
1973 [23; Nymmik, 1981 [30]; Tang, 1984 [3]]; 
Nishimura, 1999 [32J (see also earlier papers of 
the same experiment [Ml EH EE] ) • 



Figure 4. Measuremets of (e + + e )- 
spectrum, 2005-2008 and some older data: 
Nymmik, 1981 [3D], Nishimura, 1999 [32] , 
ATIC-2, 2005 [35], PPB-BETS, 2007 [Ml H], 
ATIC-2, 2008 [2J. 



of experiments containing energy points above 100 GeV. Electron spectra, measured in several 
well-known but relatively low-energy experiments such as AMS-01 [24J and HEAT |25] were 
excluded from the review. 

We divided into two classes all the available data on the high-energy cosmic-ray electron 
spectrum. We called the first class of data 'low-resolution experiments'. This class related to 
the spectral data measured with relatively large energy binning for the electron energy: more 
than approximately 12%— 15% per bin. Even if a spectrometer could provide a higher energy 
resolution, too large energy binning prevented the observation of details in the spectrum which 
were finer than the width of the energy bin. The second class of data addressed measurements 
with much more narrow energy binning in the electron spectrum together with the high energy 
resolution of the apparatus. In such 'high-resolution experiments', simultaneously, the energy 
binning and the resolution were better than 5%-8%. All available data except the only ATIC 
experiment |3j at the time of the preparation of this review, were of the first (low-resolution) class, 
and the high-resolution class of measurements was represented by the only ATIC measurements 
[3J. In this section, we review the data from the low-resolution experiments. The high- resolution 
experiment [3J is discussed separately in Section 

Figure EJ shows a compilation of the electron spectrum data obtained before 2000. The 
experiments were of a different nature: these were both balloon-borne [Ml E7J [Ml E3 Ell E2] and 
satellite [M] • They used different techniques including emulsion chambers |26[ |2"91 [30] , Cerenkov 
counters/calorimetric [2?, [28] , transition radiation/calorimetric [31j . A review of experiments 
before 2000 may be found also in [12]. The Moskalenko and Strong conventional model [IB] 
and the uniform thin disk model, described in the Section 12.31 eq. (|16p . are shown, also, in 
Figure El together with the data. The parameters of the thin disk model were 70 = 2.3, 7^ = 2.7, 
A = 1.0, Qo = 450m~ 2 sr~ 1 sec~ 1 , Q s _ e - = 24m _2 sr _1 sec _1 , Q s e + = 62 m _2 sr _1 sec _1 . Here 
Qs e+/Qs e~ ~ 2.6 was the ratio of positrons to negative electrons production rates in p-p 
collisions [16] . 

The situation at the beginning of the 2000's looked very controversial. There were two 




Figure 5. Measuremets of (e + + e )- 
spectrum: ATIC-2, 2008 [2], H.E.S.S., 2009 
[TO], Fermi/LAT, 2010 0. 



Figure 6. Measuremets of (e + + e )- 
spectrum, the current state of the art: PPB- 
BETS, 2007 [Ml P; H.E.S.S., 2009 [TO]; 
Fermi/LAT, 2010 [8J; PAMELA, 2010 0; 
ECC(Nishimura), 2011 [57]; MAGIC stereo, 
2011 mi; and ATIC2+ATIC-4 - the spectrum 
of the paper [3] shown with enlarged energy 
bins. 



different groups of data. The first group suggested a soft electron spectrum (7 > 3) at energies 
above ~ 50 GeV: Silverberg, 1973 [28]; Nymmik, 1981 [30J; Tang, 1984 [31]; and Nishimura, 
1999 [32J (see Figure [3]). The second group suggested a hard electron spectrum (7 < 3) and 
a much higher intensity of the spectrum above 100 GeV: Meyer 1971, [27] and TIFR(Anand), 
1973 [29j. The behavior of the data in the experiment Zatsepin, 1971 |26| was somewhat median 
between these two groups. The statistics of the experiments were low and it was impossible to 
discuss seriously the existence of structures like bumps in the spectra. It was difficult to obtain 
some definite conclusions on the basis of the collection of the data completed before 2000. 

There were no new data on high-energy cosmic ray electron spectrum from 2000 to 2004, but 
in 2005 the preliminary electron spectrum measured by the Antarctic long-duration balloon- 
born experiment ATIC-2 (the second flight of the ATIC spectrometer) were published [35] (see 
Figure H]) . ATIC covered the energy range from 30 GeV to more than 2 TeV and provided the 
data with much better statistics than from all previous experiments. There were a number of 
new features in the data. The electron spectrum showed that the spectral index was close to 3.0 
in the energy range 30-300 GeV; and there was a hint on energy cut-off near ITeV (however, 
with low statistical confidence) and, more importantly, a prominent bump between 300 and 
600 GeV. Thise behaviour disagreed with the conventional model; however, between 2005 to 
2007, the phenomenon was not discussed seriously. 

In 2007, the data from the Antarctic PPB-BETS experiment was presented at the 30th 
International Cosmic Ray Conference and published in 2008 [H [36] (see Figure 0]) . The PPB- 
BETS spectrum demonstrated, also, a bump-like structure between 100 and 700 GeV which was 
contradictory to the expectation of the conventional model. The statistics were lower than in 
previous ATIC experiment and the shape of the bump was somewhat different than in ATIC 
data. Near the end of 2008, the ATIC paper on the electron spectrum was published in Nature 
[2]. The bump, in the preliminary spectrum [35J, was confirmed and all this provoked a wide 
discussion of the nature of the phenomenon in the literature (see below). The bumps, observed 
by ATIC and PPB-BETS, looked similar to the bump observed, also, in 1999 in the spectrum 



measured in the Japanese emulsion experiment [32] (however, with much lower statistics, see 
Figure HJ) . Consequently, the existence of this bump- like feature in the electron spectrum proved 
rather more conclusive this time (2008). 

The situation became less clear after the publication, in 2009, of the electron spectrum 
measured by the space spectrometer Fermi/LAT [7] and then in 2010 [8] (see Figure [5]). 
Fermi/LAT provided very high statistics and confirmed ATIC data in the energy range 30- 
200 GeV by showing the spectral index to be close to 3.0. However, at higher energies, the 
Fermi/LAT spectrum did not exhibit any prominent spectral features similar to the ATIC 
bump. The Fermi/LAT instrument had a lower energy resolution than the ATIC spectrometer 
due to the thin calorimeter of Fermi/LAT against the ATIC thick for electrons calorimeter. 
However, this reason could not explain the lack of the structures in the Fermi/LAT data. 
Also, in 2009, a complete analysis (low-energy analysis together with high-energy analysis) of 
the electron spectrum, measured by a ground-based Cerenkov telescope H.E.S.S. in the energy 
range 350 GeV-5 TeV, was published in the paper [10] (the high-energy part of the spectrum 
was published previously in 2008 [38]). The H.E.S.S. telescope data confirmed, with better 
statistics, the high-energy cut-off in the electron spectrum measured by ATIC. However, due to 
the high energy threshold of the telescope, it could not confirm or disprove clearly the bump in 
the energy range 250-700GeV (Figure [5]). 

New measurements of cosmic ray electron spectrum were carried out in 2010-2012 after 
Fermi/LAT and H.E.S.S. measurements (Figure [6|). The PAMELA space spectrometer electron 
spectrum, measured in the energy range 25-640 GeV, was published in 2010 [5]. In the energy 
range 25-200 GeV, the spectrum was similar to the spectra measured by ATIC and Fermi/LAT 
and showed the spectral index to be close to 3.0. However, at energy near 200 GeV, the 
PAMELA spectrum showed some down step-like features. For energies above 200 GeV, it 
provided an absolute flux lower than the ATIC and Fermi/LAT spectra. However, the statistics 
of the data were low and this conclusion was not firm. The PAMELA collaboration published, 
also, (in two different versions), the spectrum of negative electrons [6]. However, it would be 
incorrect to compare it with the total e~ + e + spectra measured by other experiments discussed 
here. The collaboration of Japanese emulsion experiments [33} 134"! |2T| [32] published, in 2011 
[37] . the electron spectrum with updated and improved statistics. ECC was the new name 
of the experiment. This ECC spectrum confirmed well the previous version of the data [32]. 
The electron spectrum was measured, also, by the new MAGIC ground-based Cerenkov stereo 
telescope. The first preliminary version of the spectrum was published in the dissertation [39] 
and, then, improved data was published in the paper |11] (shown in Figure [6]). The spectrum 
showed a weak bump-like feature with the maximum near 250 GeV; however, the statistical 
significance of the structure was low. Finally, the electron spectrum, measured by the the ATIC 
spectrometer 4 th flight was obtained and published together with updated data of the ATIC- 
2 [3] spectrum. Generally, the ATIC-2+ATIC-4 spectrum (see Figure [(JJ) confirmed previous 
ATIC-2 [2] data with minor corrections. The paper [2] main purpose was a measurement of 
high-resolution spectrum of electrons (see Section [S]). However, for convenience and in order to 
compare it with other data, the spectrum, measured in [2], is shown in Figure [6] with enlarged 
energy bins. 

The current situation with experimental electron spectrum looks very controversial. There 
are a number of modern experiments which provide a bump-like structure in the energy range 
250-700 GeV - ECC; PPB-BETS; ATIC; and, might be, MAGIC. However, in this region, the 
Fermi/LAT spectrometer showed no bump-like feature. It is important that the Fermi/LAT 
experiment is a very high-statistical one; however, the statistical significance of the data from 
PPB-BETS and ATIC experiments also was rather good. Obviously, in some experiments, 
there were unaccounted systematic errors and, even, the expected reported systematics were 
generally large (for example, see the systematic error corridor for the ATIC-2+ATIC-4 spectrum 




Figure 7. Measuremets of positron fraction 
e + /(e + + e") before 1990: Fanselow, 1969 00]; 
Daugherty, 1975 0JJ; Buffington, 1975 02]; 
Muller & Tang, 1987 03]; Golden, 1987 03]. 
Models: Protheroe (1982) [32]; Moskalenko, 
Strong (1998) [IB]; Thin disk uniform model - 
this work, 
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Figure 8. Measuremets of positron frac- 
tion e+/(e+ + e~) before PAMELA (2008): 
Basini, 1991 05]; Golden(TS93), 1996 0B]; 
CAPRICE94, 2000 07J; HEAT-94+05+2000, 
05]; AMS01, 2007 02]. 
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Figure 9. Measuremets of positron fraction 
e + /(e + + e~): the first PAMELA data and 
some previous measurements: CAPRICE, 
2000 07]; HEAT-94+05+2000, 2004 08] ; 
AMS01, 2007 09]; PAMELA, 2009 0]. 



Figure 10. Measuremets of positron fraction 
e + /(e+ + e~), all data after 2000: HEAT- 
94+05+2000, 2004 08]; AMS01, 2007 09]; 
PAMELA, 2009 0]; PAMELA, 2010 [50] : 
Fermi/LAT, 2011 [9]. 



in Figure [6]). We can say, with confidence, that now, the main problem is not the lack of statistics 
but unaccounted systematic errors in the experiments. 

Whilst the situation in relation to the existence of some bump-like structure in the electron 
spectrum is completely unclear, the data from all high-statistical measurements are in clear 
contradiction to the conventional model prediction: the spectrum is harder than the expected 
one and has a cutoff near 1 TeV. 

3.2. The positron fraction 

Figure [7] shows the data on positron fraction in the total e + + e~ flux, obtained before 1990, 
together with the three conventional-type models. It can be seen in Figure that already the 



data of experiments in the 1970's and 1980's, showed a deviation above 10 GeV from the simple 
conventional models. There was an increase of the positron fraction instead of a reduction 
predicted by the conventional models (see eq. (|17p ). At high energies, the deviation of the data 
from the expectations was recognized clearly and, in the 1980's, noted in a number of papers 
|19l 143]. However, the statistics were too low to reach firm conclusions about it. 

Figure [8] shows the data, obtained between 1990 and 2008 (before the PAMELA experiment). 
Generally, these confirmed a deviation in the positron fraction from the conventional model 
predictions. However, the statistics were still too low and the deviation was not discussed 
widely in the literature (however, it was noted again, see [5T ] [52" | [531 154] ). 

The situation changed dramatically after the publication, in 2009 [4] (Figure E]) , of the data 
from the space spectrometer PAMELA. These high-statistical and high-precision data from the 
space magnetic spectrometer PAMELA showed definitely that the positron fraction increased at 
energies above 10 GeV. The positron fraction, measured by PAMELA as below 5 GeV, was lower 
than the data from previous experiments [471 l4"8] |4H] which provided, also, a sufficiently high 
precision in this region. This difference was not a contradiction and it could be explained by 
different charge-sign dependent solar modulations at different epochs of solar magnetic polarity. 
Such modulation variations were confirmed directly in a series of AESOP measurements [55] . 

All the modern measurements of the positron fraction after 2000 are shown in Figure [TUJ 
Recently, the PAMELA [1] 2009 data was confirmed again by the PAMELA spectrometer 
through new experimental data and with the application of three novel proton background 
models [50]. Figure flOl shows the data related to beta-fit (fit of background with using of beta- 
function, see details in [50]). Finally, the space spectrometer Fermi/LAT measured indirectly the 
positron fraction by using the Earth magnetic field to separate negative electrons from positrons 
[9]. The increasing positron fraction was confirmed up to the energies about 200 GeV, and the 
positron fraction achieved the value ~ 0.3 (see Figure fTOj) . 

Thus, there are two puzzles related to electrons in the cosmic rays: no power-law spectrum of 
electrons (possibly with a bump) (Figured]) and the increase in the fraction of positrons above 
10 GeV (Figure [TO]) . Both phenomena contradict the conventional model. Therefore, probably, 
either the conventional model is oversimplified or incorrect in some respects. The cause of these 
deviations of the data from the model should be understood. 

4. Three main ways to explain the anomalies in the electron spectrum and 
positron fraction 

The publication of the electron spectrum, measured by ATIC [2] and the positron fraction, 
measured by PAMELA [3] , resulted in a burst of theoretical thoughts which started from 2008- 
2009. Hundreds of papers with discussions and explanations of the data were published. It is 
impossible to review all the related papers here. Instead, we consider and classify here, with a 
number of examples, only three main directions of theoretical investigations. 

4-1- Conservative models 

Under the conservative way, we mean those models which try to save the conventional model 
main supposition which is to consider cosmic ray positrons as pure secondary particles. 
Generally, in such models, it is possible to explain bumps and energy cutoff in the total electron 
spectrum with sources like local SNRs or an inhomogeneity of the sources distribution [56" ] [57 ] [55]. 
However, the main problem is the PAMELA anomaly in the positron fraction. There are only 
a few papers, following the ATIC [2j and PAMELA [3] results, which try to understand them 
within the conservative model. 

The authors of paper [59J solved analytically the diffusion equation with suppositions which, 
generally, were quite similar to the ones of the conventional model \19\ 116]. However, the efforts 
concentrated on the study of the origins of the theoretical uncertainties in the estimation of 
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Figure 11. Positron fraction e + /(e + + e~): 
conservative models of T. Delahaye et al and 
N.J. Shaviv et al , and all data after 2000: 
HEAT-94+05+2000, 2004 @8]; AMS01, 2007 
US]; PAMELA, 2009 0]; PAMELA, 2010 [50]; 
Fermi/LAT, 2011 0. 
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Figure 12. A series of a single-pulsar electron 
spectra for different pulsar ages with the 
following other parametrs: r = 1 kpc, 7 = 1.3, 
E cnt = ITeV, V W = 5 • 10 49 erg, b = 
1.4 • 10~ 16 (GeV • a)" 1 , D = 3 • 10 28 cm 2 • s" 1 , 
5 = 0.3. The data of ATIC-2+ATIC-4 and 
H.E.S.S. (see Figure E|) shown, also. 



the positron flux (diffusion coefficient; cross sections; and other ISM parameters). Therefore, 
the authors of [59] obtained not only one result but, also, a number of rather wide corridors of 
predictions. In order to calculate the positron fraction, they used an experimental estimation 
for the total electron spectrum based mainly on the data of AMS-01 [24J. Figure [TT1 shows the 
corridor of predictions of [59J in the supposition of 7 = 3.35. It can be seen that the prediction 
cannot be reconciled with the data even with all the model uncertainties. It is especially true 
with the last Fermi/LAT data for the positron fraction. We would like to note that the spectral 
index of 7 = 3.35 for the electron spectrum, used in the calculations, looks overestimated as 
shown by the last data of ATIC [2], Fermi/LAT [8j, and PAMELA [5|. With a more realistic 
index 7 « 3.0, the model of paper [59] would predict an even steeper function for the positron 
fraction and the contradiction between the model and the data would be even firmer. 

In the paper [60j, it was argued that inhomogeneity of cosmic ray sources, due to the 
concentration of SNRs towards the galactic spiral arms, could explain the anomalous increase of 
the positron fraction above 10 GeV. The idea was as follows. If the observer was located between 
galactic arms, then the primary electron flux produced by SNRs, which are located mainly in the 
arms, was suppressed due to radiative cooling. This was in contrast to the secondary positrons, 
sources of which - cosmic ray protons - were distributed homogeneously in the Galaxy. The 
prediction of the model is shown in Figure [TTJ It can be seen that, actually, the model fits 
the data of PAMELA below 100 GeV. However, unknown to the authors of the paper |60j . 
it does not fit the more recent data of the Fermi/LAT observatory. We would like to note, 
also, that, actually, the Sun is not located between galactic arms. Instead, the Sun is located 
within the Orion-Cygnus arm [61]. Therefore, strictly speaking, the model is inapplicable to 
the Sun position; however, nevertheless, it points out that, generally, the sources large-scale 
inhomogeneity should be taken into account. 

Consequently, the last PAMELA and Fermi/LAT positron data seems to rule out the models 
in which the positrons are treated as purely secondary particles. Therefore, some unconventional 
source of positrons should be incorporated into the models. 



4-2. Local burst-like sources of primary electrons and positrons: pulsars and SNRs 
There is a long history concerning the idea that, in the role of local nearby sources, pulsars, 
at high energies, could determine essential features of the electron spectrum and the positron 
fraction. Starting in 1970, with C.S. Shen seminal work [13], nearby pulsars, as a possible source 
of cosmic ray electrons and structures in the electron spectrum, were discussed many times even 
before ATIC and PAMELA results (see for example [62, 52, 63J). It is important that pulsars 
produce high-energy negative electrons and positrons in equal parts (see for a review of the 
relevant physics of pulsars and pulsar wind nebulae (PWN) Appendix A of [03], and, also, 
the paper [65] and references therein). Therefore, pulsars are definitely the sources of primary 
positrons which could be relevant to explaining the PAMELA anomaly in positron fraction. 

There are, also, ideas that usual SNRs can accelerate not only negative electrons but, also, 
positrons. Suprathermal positrons, in SNRs, can be produced in beta-decay in the chains 
starting with radioactive nuclei ejected during a supernova explosion: 56 Ni, 44 Ti, 26 Al and 
22 Na [661 E3 ESI EHl E0]. A possible role of reverse shocks, in the production and acceleration 
of negative electrons and positrons, was emphasized in [691 [TO] and it was demonstrated that 
even the primary positron spectrum could be harder than the primary spectrum of negative 
electrons in [70] . It may be relevant to the PAMELA anomaly. Production and acceleration of 
secondary positrons in SNRs is considered in JTT] [72] . Independently, as regards the production 
and acceleration of positrons in SNRs, it is argued in the literature over a long time, that nearby 
SNRs could produce complicated structures in the high-energy part (~1 TeV) of the cosmic ray 
electron spectrum [56 1 157 } [58] . These may be relevant to the ATIC and PPB-BETS bump. 

In PWNe and SNRs, the character time of electrons acceleration is no more than 10-30 kyr 
and the character size of a relevant source is a few parsecs. At the same time, the age of the 
sources, relevant to our observations of electrons in cosmic rays with TeV-energies and less, was 
greater than 200 kyr and the space scale of the observations was of the order or even much 
greater than 100 pc. Therefore the approximation of point-like and instantaneous source is 
quite reasonable to estimate the electron flux from pulsars or SNRs. This approximation was 
introduced already in the C. S. Shen 1970 paper [14] . and, now, it is known as a burst-like 
approximation. 

In the burst-like approximation, the source function for a source with spectrum Q(E) located 
at r = and t = is: 

Q(r,t,E) = Q(E)5 3 (r)5(t). (18) 

By substituting the source function ()18p to the Green-function solution of the transport equation 
([9]), one obtains easily the observed density of electrons at arbitrary point r and arbitrary time 
t > in a simple analytical form: 

fE*\ 2 e -r 2 /2A 2 E 

*• *■*>-«*•>(*) ( ^xhb,e-) - e ' IE4) = t^e- < 19 > 

Here, E*(E,t) is the energy of an electron cooled down to E due to the radiative losses during 
the time t. Please note that p(r,t,E) = for E > E max (t) = l/(bot), therefore, there exists an 
exact cutoff energy E max (t) which depends on the age of the source t. 

The most natural way to understand the PAMELA anomaly and the ATIC bump is to 
consider nearby pulsars as sources of electrons since, surely, they provide large fluxes of positrons. 
To present a source spectrum Q(E) of a pulsar the exponentially-truncated power-law spectrum 
is used widely: 

Q(E) = Q E-^ exp(-E/E cnt ). (20) 

In the source function, there are three free papameters - Qo, 7, E cut - which should be selected 
on the basis of observations or from some other additional considerations. 




Figure 13. The predicted electron spectrum and positron fraction for pulsars in the ATNF 
catalog together with conventional background as calculated in [64, FIG. 4] (the figure reproduced 
by permission of the authors). 



The source spectral index 7 for a PWN, is estimated to be 1.3 =t 0.3 on the basis of radio 
astronomic data for seven pure PWNe derived from the last version of D.A. Green catalogue [73] 
(see, also, the discussion in |64j ) . It can be seen that PWN electron source spectra are harder 
than electron source spectra of SNRs (7snr ^ 2.0). 

For 7 < 2, the amplitude Qo may be estimated from the pulsar complete energy deposit to 
electrons using the following simple formula: 

POO 

/ Q(E)EdE = Q T(2 - 7 )^cut 7 = rjWo, (21) 
J 

Where Wq is the spin-down energy of the pulsar which transforms mainly to the energy of 
electromagnetic radiation and energy of electron-positron pairs, and r\ is the conversion factor 
for spin-down energy to transform to the energy of electron-positron pairs. 

There are significant uncertainties in the estimations of Wo for known pulsars. However, 
generally, it is assumed that, for the majority of pulsars, Wo is between 5 • 10 48 and 5 • 10 50 erg 
[M] . [741 Table 4]. The conversion coefficient 77 is expected to be ~ 1 but, effectively, it could 
be somewhat less (down to ~ 0.1) due to energy losses during the retention of electrons within 
PWN before destruction of the PWN and ejection of the electrons into the ISM [64j . Therefore 
10 48 -10 50 erg is a reasonable estimate for tjWq. 

The cutoff energies E cnt are not known well. However, there are indications that they might 
be as high as several tens of TeV. For example, the electron cutoff energy for Vela PWN, as 
estimated by the gamma-telescope H.E.S.S., is 67 TeV [75J 4 ! 

By variation of the parameters r, t, Qo,^y, E cut , one could obtain a variety of different single- 
pulsar spectra with various shapes and intensities. Figure [12] shows an example of a series of 
a single-pulsar electron spectra with different ages but with the same other parameters. It can 
be understood from Figure [12] that within the existing uncertainties of parameters of known 
nearby pulsars the intensity and shape of the spectrum provided by a combination of several 
pulsars together with conventional background can be agreed with the data of ATIC, Fermi/LAT, 
H.E.S.S. etc. by a variety of ways. Really, many recent papers describe the data by this manner 
(see for example [73 [771 IM1 E3 EH [78] ; for review see [79]). 

4 However, the estimated spectral index from gamma-ray flux is 7 ~ 2.0 and Vela is not pure PWN according to 
the catalogue [73], it is PWN+SNR nebula. 



In order to describe the electron spectrum high energy features in such models, an essential 
part or, even, the most of the flux, at energies E > 200 GeV, is given by a contribution of nearby 
pulsars with equal parts of positrons and negative electrons. Therefore, generally, the positron 
fraction increases from the low values like ~ 0.05 at E ~ 10 GeV up to the values > 0.2 at 
energies 300-500 GeV. Consequently, such models explain simultaneously the PAMELA positron 
anomaly together, very naturally, with the no-power-law behaviour of the electron spectrum. 
This result was demonstrated in all cited above papers [TBI E7J EH E3 EH EE] and, also, in the 
earlier works |51l [52] . The calculations for all pulsars in the ATNF catalogue carried out in the 
Malyshev et al. paper [M], are shown as an example in Figure [131 

The possibility of SNRs to accelerate not only electrons but, also, positrons does not interfere 
with working models of pulsars as sources of electrons. Moreover, for the same reasons as pulsar 
models, the acceleration of positrons only by SNRs, without pulsars at all, can explain, in 
principle, the electron spectrum and the positron fraction simultaneously [71j . In all these 
cases, the models kernel part is the possibility of local instantaneous sources to accelerate both 
negative electrons and positrons. 

4-3. Dark matter annihilation and decay 

WIMPs annihilation (or decay) is the third main way in understanding the features of the 
electron spectrum and the positron fraction anomaly. WIMP means 'weakly interacting massive 
particle' - one of the most promising candidates for the DM (see for review [80} \81\ [821 183j). 
For the first time the possible anomaly in the positron fraction at energies above 10 GeV was 
directly attributed to the annihilation of WIMPs in the paper of A. J. Tylka [M] of 1989. Later, 
even before 2000, this question was studied several times [851 [86, 87, 88J and hundreds of papers 
were published following ATIC and PAMELA results. In the electron spectrum, the bump-like 
feature was considered to be a possible signature of WIMPs annihilation to the electron-positron 
pairs already in the PPB-BETS paper [lj. However, it was noticed that the statistics were too 
low to make definite conclusions. For the first time in the electron spectrum, the observed bump 
was attributed clearly directly to a possible signature of the annihilation of DM particles in the 
ATIC paper [2]. The mass of the particle was estimated as ~ 620 GeVH. Now there are numerous 
papers which consider the no-power-like shape of the electron spectrum by similar ways (see for 
review [79] ) 

Many-body simulations show that WIMP-like DM halo of the Galaxy must be clumpy |89j . 
In the paper [90J, it was noted, for the first time, that the DM clumps (subhalos) could enhance 
many times the DM annihilation rate, and, therefore, the clump of the DM was important for 
the problem of cosmic ray electrons ejection. Consequently, there are two different but connected 
problems: firstly, the observation of electrons produced by the smooth Galactic DM halo and, 
secondly, by nearby local clumps. We considered two distinct analytical solutions of the electron 
transport equation ([8]) related to these problems. 

4-3.1. Smooth Galactic DM halo. In order to obtain a simple analytical solution related to a 
smooth Galactic halo, we supposed the source of electrons to be infinitely large in space and 
homogeneous, with the source spectrum Q(E): 



By substituting the source function (|22p to the general Green function solution (|9j) of the 
transport equation, one obtains easily the solution for the observable electron spectrum: 



Q(r,t,E) = Q(E). 



(22) 




(23) 



5 It was done in [2] for the case of Kaluza-Klein particles which produce one electron-positron pair per annihilation. 




Energy, GeV 



Figure 15. Same as Figure Q3]but with boost 
factor 500. 



Figure 14. A typical prediction of a 
simple analytical model (|25p for the electron 
spectrum produced by homogeneous DM-halo 
without boost factor. Dashed line - the 
spectrum of electrons from DM annihilation, 
the solid line - the sum of DM spectrum 
and conventional background. Dark matter 
peak is invisible above the background. The 
parameters of the model are: Mwimp = 
1.1 TeV, k = 4, (av) = 3 • 10 _26 cm 3 s _1 , 
n = 0.3 GeV • cm _3 /M W i M p, Qo = 4.5 • 
10~ 33 cm~ 3 s -1 . 
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Figure 16. The predicted positron fraction 
for the Figure [TBI model. 



The energy spectrum of electrons, ejected in the annihilation of WIMPs, could be delta-like in 
the simplest case or, in general, more complicated. We considered only the simplest delta-like 
case. The results for more general situations could be understood easily on the basis of the 
simplest one. The source spectrum for delta-like annihilation spectrum is 

Q{E) = Q S(E - E ), (24) 

and from (|23p one immediately obtains 

f ^- E<E Q 
p(E)={ b E^ tj<tj * (25) 

[ 0, E > E 

The rate of annihilation Qo m the formulae f)24[) and (|25p may be estimated as 

Qo = hn 2 (av), (26) 

where k is the number of electrons per one annihilation; n [cm -3 ] is the WIMPs density in 
the smooth Galactic halo; a is the cross section of annihilation; and v is relative velocity of 
WIMPs. Here, we assume, also, for simplicity that the WIMP particle is its own antiparticle 
(like Majorana fermions). The density n is estimated as n ~ 0.3 GeV • chi~ 3 /M\vimp; and, 



for the WIMPs annihilation cross section, it is known from the cosmological abundance of DM 
(av) ~ 3 • 10 _26 cm 3 s _1 |81j . This value is expected to be independent on v for annihilation in 
S- mode (the Bethe law). 

As an example, we consider here the annihilation of WIMPs to two electron-positron pairs 
and suppose Mwimp = 1-1 TeV. Then k = 4 and n = 2.7 • 10~ 4 cm~ 3 . From (j26|) one obtains 
easily Qq = 4.5 • 10 -33 cm _3 s _1 and the respected spectrum calculated with equation (|25p is 
shown in Figure UM 

It can be seen that the amplitude of the DM spectrum is almost three orders of magnitude 
less than the amplitude of experimental electron spectrum and that the DM deposit is even 
invisible above the conventional background. This is due to too low cross section value 
(av) ~ 3 • 10 _26 cm 3 s _1 defined by the cosmological data and it is generic feature of all DM 
models of an origin of cosmic ray electrons. In order to obtain reasonable amplitude for the 
DM annihilation deposit to the electron spectrum, one should suppose some mechanism of 
increasing the annihilation cross section (av) at low energies - the so-called cross section boost 
factor. Possible mechanisms may be either a Sommerfeld effect or a Breit-Wigner enhancement 
factor (see for review [91]). Figure [T5l shows the DM spectrum with the same parameters, as 
described above (see also Figure fbl|) but with a 500 boost factor. In Figure \15\ the very sharp 
peak is not a problem of this model but only an artifact of a delta-like annihilation spectrum 
of electrons (|24p . For a more smooth annihilation spectrum, the observed spectrum would be a 
simple convolution of the peaks like in Figure [15] with the actual source annihilation spectrum 
and the shape of the predicted spectrum being smoother. 

Negative electrons and positrons are produced in the annihilation of WIMPs in equal parts. 
Therefore, if some DM model predicts no-power-law spectrum or bumps in high-energy part 
of the electron spectrum with sufficiently high amplitude such as in Figure [151 it produces 
almost automatically reasonable behaviour of the positron fraction, as shown in Figure [161 This 
important feature of the DM models does not depend on details like the annihilation spectrum 
etc., and the feature is very general. The origin of this feature is similar to the origin of the 
similar feature of the pulsars models (see section 14.21 and Figure fT3j) . Many papers ( [92], [931 EI] 
and others, for review see [79]) demonstrated the possibility of smooth halo DM models to 
predict, at the same time, bumps in the electron spectrum and the reasonable positron fraction 
consistent with the PAMELA anomaly. 

4-3.2. Subhalos (clumps) From the late 1990's, it was known from numerical many-body 
simulations [9i\ [95] that the Galactic dark mater halo had to be clumpy since there was an 
incomplete merging of primordial DM substructures which formed the Galactic halo. Recent 
many-body simulations, with the code Via Lactea II [89], showed that the subhalos (clumps) 
were distributed with approximately equal total mass per decade of subhalo masses over the 
range 1O 6 M0-1O 9 M0, where M & was the Solar mass. Within the sphere with a 400 kpc radius 
with the center at the core of the Milky Way, it is expected that there is ~ 10 subhalos of 
~ 1O 9 M0, ~ 10 4 subhalos of ~ 10 6 Mq etc. The mean distance, between 10 6 Mq subhalos, is 
~ 10 kpc and, therefore, it is approximately the expected distance from the Sun to the nearest 
subhalo of this mass. The mean DM density, within such a clump, is expected to be hundreds 
of times greater than the mean density of the smooth Galactic halo, and the density near the 
central parts of a clump (a cusp) - is expected to be thousands of times greater. The annihilation 
rate is proportional to the square of the DM density. Therefore, there is a natural 'boost factor' 
for the annihilation rate and the related flux of electrons may be rather high near the Sun. This 
problem should be investigated carefully. 

The scale of sizes of ~ 1O 6 M0 clumps is ~ 100 pc, the scale of sizes of most dense and bright 
cusps of clumps is only ~ 10 pc [89 1 . A point-like approximation is quite reasonable for the 
observations of the clump from distances ~ 1 kpc and greater. Also, in time, the related source 
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Figure 17. A family of electron spectra 
produced by a single DM clump without boost 
factor for different distances from an observer 
to the clump. The total mass of the clump is 
about 2.5 • 1O 6 M0, the total annihilation rate 
is Qo = 1.7 • 10 34 s _1 with delta-like spectrum 
of electrons concentrated at 550 GeV. 



of electrons should be considered as permanent and constant. Therefore, supposing as above a 
delta-like spectrum for annihilation electrons, we consider the source function to be: 

Q(r, t, E) = Q 5 3 (r)5(E - E ). (27) 

Substituting the source function ([2?]) to the Green function solution (J9]) of the transport equation 
one obtains easily a simple analytical solution which helps to understand the DM clumps physics 
related to cosmic ray electrons: 

exp[-r 2 /2A 2 (e,£ )] 
p{r > t) = Q ° b EH2^\3(E,E y (28) 

In order to estimate an observed electron spectrum from a dark-matter clump, we supposed 
a simple clump toy model. The clump was represented by two concentric spheres. The outer 
sphere presented a volume of the clump; the inner sphere presented a cusp. For the volume, 
we supposed pdm (volume) = 2.5M /pc 3 , ^(volume) = (lOOpc) 3 ; for the cusp we supposed 
Pdm(cusp) = 45M /pc 3 , V(cusp) = (10pc) 3 . The total mass of such clump was about 
2.5 • 1O 6 M0. These parameters looked possible for a clump with such sizes (see [E2]); however, 
the density and mass could be, also, a ~ 5 times less [89J). Therefore, our toy model clump was 
relatively dense and bright and our estimates for the amplitude of the electron spectrum should 
be considered as very optimistic. 

For this toy model, with the parameters of DM particles as in the Section 14.3.11 one obtains 
easily Qo = 1.7 • lO^s" 1 in the eq. (|27p . Figure [TTI shows a family of the spectra for various 
distances from the Sun to the clump. It can be seen that, even without a boost factor, the 
deposit from the nearby DM clump to the electron spectrum may be observable, in principle, 
but only for a very close clump r < 0.1 kpc (in fact - within the clump). This situation is very 
unrealistic. For more realistic distances r ~ 1 4- 10 kpc the clump deposit is negligibly small 
compared to the observable flux of electrons and it may be observed only in the case of high 
boost factors like 10 2 4- 10 3 . 

If, actually, the annihilation energy spectrum of the electrons is delta-like, then, it is possible, 
in principle, to obtain a very sharp peak in the electron spectrum. This situation may take place 
in the case of a DM clump located very close to the observer (see Figure fT7|) . However, we would 
like to emphasize that it is impossible, in principle, to have two or more sharp peaks related to 
DM at different energies and at the same time. As clear from Figure [TTI if there are two or more 
clumps at different distances from the Sun, they must produce peaks at different energies (due 
to cooling of the electrons). However, only the nearest clump could produce a sharp peak, the 
others produce only very smooth distribution. It was noted in [6l] and direct calculations, for 
many-clump structures, revealed, also, this feature [961 197]. 
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Figure 18. Toatal electron flux and positron fraction from boosted annihilating DM in the 
host halo plus 1 kpc-clump as calculated in |98} FIG. 7] (the figure reproduced by permission of 
the authors). See details of the calculation in |98j. 



A possibility to predict simultaneously bump-like structures in the electron spectrum and 
rising positron fraction, is a generic property of subhalo DM models. This is an implication of 
the equal fluxes of negative electrons and positrons in the annihilation of DM particles and it 
is quite similar to the same possibility of the pulsars models and DM smooth halo models (see 
above). However, the spectrum, produced by a DM clump, is generally harder than the spectrum 
produced by a smooth halo. The result is a deficit of positrons near the energies ~ 100 GeV 
relative to PAMELA data for a pure DM clump scenario [UB]. However, the dark matter smooth 
halo, together with a deposit of clumps, can describe bumps in the electron spectrum and, at the 
same time, the anomaly in the positron fraction. It was demonstrated in many papers (see for 
example [Ml Ell SB] , for review, please see [Z2]). As an example of such predictions, Figure [TBI 
shows the total flux of the host smooth halo, together with a nearby 1 kpc clump calculated in 
the paper [98] . 

We did not discuss in detail the possibility of explaining the ATIC and PAMELA anomaly 
by decay of DM particles with generation of electron-positron pairs. Such a possibility is 
discussed widely in the literature (see for example [99 J and references therein) . This explanation 
shares many features of DM annihilation with the exception of a strong boosting due to a high 
concentration of DM in local clumps. Particularly, DM decay cannot produce several distinct 
sharp peaks in the electron spectrum as well as DM annihilation cannot, as was discussed above^l 

4-4- To the model selection 

It can be seen from the review, given in Section 0] that the conservative way could not explain the 
positron anomaly; therefore, it was ruled out. On the contrary, both the pulsars scenario and the 
DM scenario could explain successfully the data of low-resolution experiments (see definition of 
low-resolution experiments in Section 3.1). Therefore, it was impossible to select between these 
two explanations only on the basis of low-resolution electron data. 

In the DM scenario for electrons, there are several troubles. A natural prediction of DM 
models is an anomaly in the antiproton fraction in cosmic rays similar to the positron anomaly; 
however, a similar anomaly was not observed in the PAMELA experiment |1UU|, HOlj . This 



In order to produce several distinct sharp peaks in the spectrum of electrons, the only possibility is to have, 
at the same time, 1) DM decay or annihilation with several distinct delta-like peaks in the energy spectrum of 
electrons; and 2) location of the Sun inside a DM-clump, close to its cusp. It looks absolutely improbable. 




Figure 19. The ATIC spectrometer 
(ATIC-2 configuration). 1-silicon matrix, 2- 
scintillator hodoscopes, 3-carbon target, 4- 
BGO calorimeter. 



Figure 20. Total BGO energy deposit for 
150 GeV (left) and 300 GeV (right) electrons 
obtained in the beam tests [106 J . 



trouble is avoidable by suggestion of a leptophilic annihilation of DM particles |102j . The 
strongest constrains on the DM annihilation models come from gamma-ray astronomy. There 
are large uncertainties associated with model-dependent or poorly known astrophysical factors 
like DM density profiles. Nevertheless, the obtained upper limit to {av) excludes the smooth 
Galactic halo as the only source of the ATIC and PAMELA anomaly [103] 1104"] . However local 
DM substructures like clumps are not excluded by gamma-ray constraints since they may in 
principle produce strong signal even without boost factor (see Figure [TT|) . 

However, local DM substructures, like clumps, are not excluded by gamma-ray constraints 
since, in principle, they may produce strong signals even without a boost factor (please see 
Figure fT7|) . Therefore, DM models cannot be excluded and the nature of electron anomalies 
remains unclear. 

The situation may be changed dramatically if the results of high-resolution ATIC experiment 
[3], considered in the next section, are confirmed. 

5. The fine structure of the electron spectrum measured by the ATIC 
spectrometer 

The ATIC (Advanced Thin Ionization Calorimeter) balloon-borne spectrometer was designed 
to measure the energy spectra of nuclei from hydrogen to iron with an individual resolution of 
charges in primary cosmic rays for energy region from 50 GeV to 100 TeV. It was shown that 
ATIC was capable of not only measuring the spectra of cosmic ray nuclear components but, 
also, the spectrum of cosmic ray electrons [105j . In order to separate the electrons from a much 
higher background of protons and other nuclei, the differences in the shower development in the 
apparatus for electrons and nuclei are used. The low-resolution spectrum (see the definition of 
the term in Section I3.1j) , measured by ATIC in this way, was published in the paper [2] . 

The ATIC apparatus is comprised of a fully active bismuth germanate (BGO) calorimeter; 
a carbon target with embedded scintillator hodoscopes; and a silicon matrix that is used as a 
main charge detector as it shown in Figure [191 The ATIC project had three successful flights 
around the South Pole in 2000-2001 (ATIC-1), 2002-2003 (ATIC-2) and 2007-2008 (ATIC-4). 
ATIC-1 was a test flight and is not discussed here. The details of construction of the apparatus 
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Figure 21. The spectrum of electrons with- 
out the subtraction of proton background and 
without an atmospheric correction as mea- 
sured in the ATIC-2 and ATIC-4 experiments 
[3J. The size of energy bin is 0.035 of the dec- 
imal logarithm of energy (8.4%). 
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Figure 22. Total statistics of the electron 
spectrum in ATIC-2+ ATIC-4 flights in each 
energy bin. In this Figure, the summary 
spectrum corresponds to the spectra in 
Figure ED In the range of 200-810GeV, 
there were 701 events studied for statistical 
significance of the structure. 



and the procedures of its calibration are given in [1071 11081 \Wb\ [TPl?] . 

The calorimeter of ATIC was thick for electrons. It was 18Xq [Xq is a radiative unit) in 
ATIC-2 and 22Ao in ATIC-4. As a result, the energy resolution for electrons of the instrument 
was high in the total energy range of 30 GeV-1 TeV. The simulation predicted and the beam 
tests confirmed that the resolution, in terms of half width on the half height, was not worse than 
3% |106j (see Figure [20]) . Such high resolution, together with relatively high statistics, allows 
a meaningful measurement of the electron spectrum in the high-resolution mode - with energy 
bins as small as 5-8% (0.020-0.035 of decimal logarithm of energy). 

Such high-resolution measurements were carried out for ATIC-2 and ATIC-4 flights in [3] and 
revealed a fine structure in the region of "ATIC bump" (200-600 GeV) which was reproduced in 
both flights, see Figure ED It can be seen that there were three peaks near 250 GeV, 350 GeV 
and 450 GeV respectively and three dips near 220 GeV, 300 GeV and 420 GeV. Figure EH 
shows the spectra of electrons without subtraction of the residual proton background and 
without an atmospheric correction. Both the proton background and the factor of atmospheric 
correction are smooth functions and can not disturb the fine structure notably. Only these 
'raw' spectra are the input for the estimation of the statistical significance of the fine structure. 
The statistical significance was estimated in [3] by two different ways. The first way was to 
use the usual x 2 -test for the total ATIC-2+ATIC-4 spectrum. The second one was to use the 
degree of correlation between the structures measured separately by ATIC-2 and ATIC-4. In 
both ways the statistical significance of the fine structure was estimated as 99.7%. Please note 
that the bin statistics are not really low in the spectrum - it varies from 70 events/bin to 40 
events/bin near the three mentioned above peaks of the fine structure for the total spectrum 
ATIC-2+ATIC-4, see Figure E2J The observed fine structure is expected to be stable against 
all known possible systematic errors since they cannot produce any short-scale structures in 
the spectrum. Moreover, the observed fine structure passed and survived a number of tests on 
stability against possible systematics [3J. 

Although the statistical significance of the observed phenomenon is not low and it looks 
systematically firm, it should be confirmed by independent experiments to be considered 



absolutely seriously of course. However, if the structure is confirmed, it is very important 
in understanding the nature of the electron spectrum anomalies and the nature of DM. As we 
explained in Section 14.3.21 above, the annihilation of DM particles could not produce several 
narrow peaks in the electron spectrum. Therefore, DM annihilation could not be an origin of 
the observed fine structure. Moreover, since the amplitude of the fine structure is very high, 
most of the intensity of the observed electron spectrum, at energies above 200 GeV, could not 
be attributed to DM annihilation or decay, therefore a high boost factor for DM annihilation is 
excluded. 

The most natural explanation of the fine structure is a contribution of several nearby pulsars 
to the cosmic ray electron flux. In fact, a fine structure of the electron spectrum which was very 
similar to the observed in the ATIC experiment, was predicted on the basis of nearby pulsars 
model in the article [B3] independently on the actual observation. The predicted fine structure 
is seen clearly in Figure [T3j Moreover, it was argued in [M] that such features in the electron 
spectrum at high energies, would suggest strongly a pulsar origin of the anomalous contribution 
to the electron flux against its dark matter origin. Let us note that the authors of the ATIC 
paper [3] observed the fine structure also independently of the paper [63]. The first message 
about the fine structure was published in |110j and the first preliminary spectra were shown in 
the paper [lllj without references to [64] . 

6. Fitting the fine structure by several single-pulsar spectra 

It can be seen from Figures [21] and [22] that the observed fine structure of the electron spectrum 
shows an even higher amplitude than the structure predicted in [63] (see Figure fl~3|) . This poses 
a question: is it actually possible for nearby pulsars to produce a structure with so high an 
amplitude? This question was discussed, for the first time, in the paper [112] . There, it was 
argued that it was possible, in principle; however, an unusual supposition had to be done about 
the cutoff of the conventional electron background near 200 GeV. Now, we would like to discuss 
the same question, but supposing the usual power-low behaviour of the conventional background 
up to the energies of a few TeV. 

In order to obtain a structure of the electron spectrum, like that observed in the ATIC 
experiment, obviously, the single-pulsar peaks, in the energy range 200-700 GeV, must be very 
sharp and narrow. What may a single-pulsar spectrum look like? Five series of single-pulsar 
electron spectra, for various parameters of pulsars and ISM, were obtained with the expression 
for a single-pulsar spectrum ([19]) . These are shown in Figures [121 E31 G31 ES and [26] It can 
be seen from these Figures, that sharp and sufficiently strong single-pulsar electron peaks can 
be obtained from a variety ways with quite natural parameters. By combining these ways, the 
measured fine structure can be fitted, also, to many different ways. 

Our task was not to identify each observed feature in the measured fine structure of the 
electron spectrum with some definite known pulsar. Our opinion was that this problem was too 
difficult now for two reasons. Firstly, the distances to the nearby pulsars were known poorly 
(some other parameters, too, were known poorly). Secondly, since the birth velocities of pulsars, 
due to asymmetric supernova explosions, were very high, the pulsars known locations had little 
in common with those locations where the pulsars had ejected the electrons. The mean birth 
velocity was 450 ±90 km/sec and a velocity above 1000 km/sec was not rare [113] . Typically, the 
difference between a current pulsar location and a location of burst-like ejection of the electrons 
by this pulsar may be hundreds of parsecs, and it is exactly the scale of distances from nearby 
pulsars to the Sun. Therefore, the distances to actual bust-like sources are almost completely 
unknown and it is meaningless to try to fit the current positions of pulsars. Consequently, we 
tried to demonstrate only that the fit of the fine structure by pulsars was possible with some 
reasonable suppositions. 

Figure [27] shows an example of a fit of the electron spectrum generated by several pulsars 




Figure 23. A series of a single-pulsar 
electron spectra for different E cnt with other 
parameters r = 250 pc, 7 = 1.3, rjWo = 
1 • 10 50 erg, b = 1.4 • 10" 16 (GeV • s)" 1 , D = 
3- 10 28 cm 2 -s -1 , 5 = 0.3. 



r = 50 pc 

r = 100 pc 




Energy, GeV 



Figure 24. A series of a single-pulsar electron 
spectra for different source index 7 with other 
parameters r = 250 pc, E cn t = 10 TeV, tjWq = 
1 • 10 50 erg, 6 = 1.4 • 10~ 16 (GeV • s)" 1 , D = 
3-10 28 cm'-s" 1 , 6 = 0.3. 

— 8 = 0.3 




— 8 = 0.4 

— 8 = 0.5 

— 8 = 0.6 
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Figure 25. A series of a single-pulsar electron 
spectra for different r with other parameters 
7 = 1.3, E cut = 10 TeV, r,W = l- 10 50 erg, 
6 = 1.4 • 10~ 16 (GeV • s)~\ D = 3 • 10 28 cm 2 • 
s" 1 , 5 = 0.3. 



Figure 26. A series of a single-pulsar electron 
spectra for different 5 with other parameters 
r = 250 pc, 7 = 1.3, ^ut = 10 TeV, rjW = 
1 • 10 50 erg, 6 = 1.4 • 10~ 16 (GeV • s)" 1 , D = 
3 • 10 28 cm 2 • s _1 . 



with the conventional background. The conventional background is exactly the same as in 
figures [3HS The propagation parameters for this picture are 5 = 0.3, Dq = 3 • 10 28 sm 2 s _1 , 
60 = 1.4- 10 _16 (GeV-s) _1 . The points in the figure denote the high-resolution electron spectrum 
measured by ATIC-2+ATIC-4 with the correct absolute normalization. The proton background 
is subtracted; and the scattering of electrons in the residual atmosphere above the apparatus 
is corrected (see [3] for details). The energy bin size is 0.035 of the decimal logarithm (8.4%). 
Below 200 GeV, the data are fitted by two distinct pulsars; however, it should be considered only 
as an illustration. Also, this region could be fitted by a continuum of far pulsars, as was done 
in paper [64J. The generated electron peaks is very sharp within the simple analytical model 
(|19j) . Actually, the single-pulsar peaks are smeared out of course by the ISM inhomogeneity. 
However, this smearing is expected to be about only 5% for E ~ 1 TeV electrons [63]. Such a 
small smearing could not change notably the obtained picture. 




Energy, GeV 



Figure 27. An example of a fit of 
the measured fine structure in the electron 
spectrum by several nearby pulsars with 
conventional background. The papameters of 
the pulsars are shown in the Table [H the 
propagation parameters are 6 = 0.3, Dq = 
3 • 10 28 smV 1 , b = 1.4 • l(T 16 (GeV • s) _1 . 
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Figure 28. The positron fraction calculated 
for the model of Figure [271 The experimental 
data are: PAMELA, 2010 [50]; Fermi/LAT, 
2011 [0J. 



Table 1. Parameters of pulsars used to generate Figure [271 
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There are many other reasonable sets of parameters which, equally well, could fit the data - 
this fit is a typical ill-defined problem. However, it is important that, in any successful fit of the 
data, the fraction of the pulsars electron flux in the total flux of electrons above 200GeV, must 
be very high to generate a structure with sufficiently high amplitude. It suggests immediately 
that the positron fraction must rise very rapidly along the energy and must reach values as 
high as ~ 0.4 near 400-500 GeV. It is illustrated in Figure [28] that was generated by the same 
sets of parameters as Figure [271 It is a very important and generic prediction of the model; 
we formulated it already in paper [112] . This prediction was done before the positron fraction, 
measured by Fermi/LAT [9], was published, and, now, the model prediction was confirmed by 
the data up to the energies about 200GeV. 

Anothr important prediction of the model is a fine structure that can be seen in the positron 
fraction curve at highest energies in Figure [28] It is not an artifact of the model. It is a generic 
implication of a variation of the pulsars flux in relation to the smooth conventional background 
in the energy range of the observed fine structure 200-600 GeV. In principle, this fine structure 
may be measured, also. However, since the amplitude of the structure in the positron fraction 
is less that the fine structure amplitude in the electron spectrum, it is a more difficult problem 
than the observation of the fine structure in the total electron spectrum. Also, there is a need 



for high-resolution measurements of the positron fraction to observe the phenomenon. 
7. Conclusions 

As understood from this review, there was a great progress, over the last ten years, in 
measurements and understanding of electrons and positrons in cosmic rays. However, there 
is an extremely urgent need for new precise and high-statistical experiments. 

The information from recent measurements of the total electron spectrum is very 
contradictory. There are a group of experiments which show a bump-like structure in the 
region 200-700 GeV: PPB-BETS PQ, ATIC [2], ECC [37], may be, also, the preliminary data 
of MAGIC Definitely, the Fermi/LAT spectrometer does not confirm the bump [8]. The 
result of H.E.S.S. telescope is unclear in this respect, due to high energy threshold [10], and 
the PAMELA spectrum [5] has too low statistics for definite conclusions in this region. Also, 
the absolute flux differs for different modern experiments on a factor up to 3 in the region 
300-500 GeV. 

The positron fraction was measured by PAMELA [U [6] and Fermi/LAT |114| 1115] up to 
200 GeV only. Moreover, the data in the region 100-200 GeV are indirect since they were 
obtained with the help of the Earth Magnetic field, and not directly with a spectrometer. There 
is an urgent need for high precision measurements of positron fraction for energies up to 500 GeV. 

The fine structure of the electron spectrum, measured by high-resolution ATIC experiment 
PJ, could not be compared with any existing data since the ATIC spectrum was the only up 
to date high-resolution measurement of the electron spectrum. The fine structure, measured by 
ATIC above 200 GeV, could not be resolved by using energy bins as wide as in the Fermi/LAT 
spectrum [8] and in all other published data. Moreover, the ATIC spectrum relates to the 
part of the Southern sky with the declination between —45° and —90°, whilst the Fermi/LAT 
spectrum integrates the flux over the whole sky. Taking into account the possible anisotropy of 
the electron energy spectrum, it would be more correct to compare the ATIC spectrum with the 
spectrum measured by some high-resolution and high-statistical experiment for the same part 
of the sky. 

The above noted lack of data is especially dramatic since a fine structure in the electron 
spectrum (and, also, in the positron fraction) may be a critical test in selecting dark matter or 
local sources like nearby pulsars as an origin of anomalies in the electron spectrum and positron 
fraction (see Section [5]). Usual experiments, with large energy binning in the spectrum, could 
not distinguish between dark matter and nearby pulsars as sources of cosmic ray electrons. 
Therefore, there is a need not only for high-statistical and safe against systematic errors 
measurements but, also, for high-resolution measurements - measurements with narrow energy 
bins ( «8% or less). 

The author is grateful for discussion with D. Malyshev, N.V. Sokolskaya, V.I. Zatsepin. The 
work was supported by RFBR grant number 11-02-00275. 
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